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Abstract 

The  paper  presents  a new  dynamic  SGS  model  of  two-way  coupling  for  large 
eddy  simulation  of  particle-laden  turbulent  flow.  The  advantage  of  this  new 
model  is  that  coupling  of  fluid-particles  SGS  components  was  taken  into  account 
and  at  the  same  time  the  coefficient  of  proposed  SGS  model  can  be  optimized  by 
Germano’s  (1991)  dynamic  procedure0'.  To  investigate  the  capability  of  this 
model,  numerical  simulations  of  particle-laden  turbulent  flow  at  Re=644  in  a 
vertical  channel  were  performed  using  this  new  dynamic  SGS  model.  By 
comparing  the  calculation  results  with  that  using  single-phase  SGS  models  which 
didn'  t consider  coupling  of  SGS  components  the  role  of  particles  SGS 
components  played  in  the  turbulence  modulation  of  fluid  flow  was  clarified  for 
particle-laden  channel  turbulent  flows.  In  addition,  the  subgrid-scale  stresses 
obtained  using  the  proposed  model  vanish  at  wall  boundary,  and  have  the  correct 
asymptotic  behavior  in  the  near-wall  region  of  the  turbulent  boundary. 

1.  Introduction 

In  large  eddy  simulation  of  single-phase  turbulent  flows,  one  major  deficiency  of 
the  Smagrinsky'2'  subgrid-scale  stress  models  is  their  inability  to  represent 
correctly  with  a single  universal  constant  to  different  turbulent  fields  in  rotating 
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or  sheared  flows,  near  solid  walls,  or  in  transitional  regimes.  Usually  additional 
modifications  to  the  Smagorinsky  model  were  made  in  the  near-wall  region  to 
force  the  subgrid-scale  stresses  to  vanish  at  the  wall  boundary  with  a standard 
Van  Driest(3)  damping  function.  The  dynamic  SGS  model  proposed  by 
Germano(1)  et  al.  (1991)  overcomes  these  deficiencies  by  locally  calculating  the 
eddy  viscosity  coefficient  to  reflect  closely  the  sate  of  the  flow.  In  large  eddy 
simulation  of  multi-phase  turbulent  flows,  it  is  being  expected  that  the  fluid 
turbulence  modulation  SGS  model  can  be  evaluated  by  Germano’s  dynamic 
procedure  as  in  the  single-phase  turbulent  flows. 

In  particle-laden  turbulent  flow,  a kind  of  multi-phase  turbulent  flows  that 
occur  in  a wide  range  of  engineering  and  scientific  research,  interaction  of 
particles  and  gas-phase  turbulent  carrier  flow  is  a problematic  research  topic  of 
both  fundamental  importance  and  practical  interest.  In  addition,  it  is  the  most 
interesting  problem  in  developing  numerical  simulation  models.  A two  way 
coupling  SGS  model  was  present  by  Yuu(4)  (1997),  in  which  coupling  of 
fluid-particles  SGS  component  was  considered.  However,  the  model  coefficient 
was  simply  decided  in  the  same  way  as  previous  works  for  the  single-phase 
turbulent  flows. 

We  present  here  a new  dynamic  SGS  model  for  two  way  coupling  LES  of 
particle-laden  turbulent  flow  based  on  Yuu'  s SGS  model.  The  couplings  of 
fluid-particles  SGS  component  on  fluid  turbulence  modulation  in  LES  was  taken 
into  account  and  at  the  same  time  the  coefficient  of  proposed  SGS  model  can  be 
dynamically  decided  by  assuming  that  the  small  scales  are  in  equilibrium,  so  that 
energy  production,  dissipation  and  interaction  of  fluid-particles  are  in  balance. 

The  main  objective  of  this  paper  is  to  introduce  this  model  and  at  the  same 
time  to  report  the  results  of  investigation  about  the  capability  and  limitation  of 
present  SGS  mode  in  predicting  the  fluid  turbulence  modulation  for 
particle-laden  turbulent  flows  through  numerical  simulations,  which  were 
performed  in  downward  particle-laden  turbulent  flows  at  Re=644  in  a vertical 
channel  using  Van  Driest  wall  function  model  and  the  proposed  new  dynamic 
SGS  model. 


2.  Improvement  of  Two  Way  Coupling  SGS  model 

For  dilute  particle-laden  turbulent  flow,  in  which  particle  volume  fraction  is 
very  small  but  the  particle  mass  loading  is  large,  momentum  equations  governing 
transport  of  the  large  eddies  was  oriented  by  filtering  the  incompressible 
Navier-Stokes  equations  with  the  particle  source  term,  where  the  Leonard  and 
Cross  terms  were  neglected. 
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The  effect  of  the  SGS  on  the  resolved  eddies  in  EQ  (1)  is  represented  by  SGS 
stress  tv  = UjUj  -ufij . In  this  work,  tv  is  parameterized  using  an  eddy  viscosity 
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hypothesis  expressed  by  EQ  (2),  and  the  particles  turbulent  dispersion  flux  terms 
n ir  and  n upj  adopt  gradient  dispersion  model  as  presented  in  EQ  (3)  and  EQ 
(4). 
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Here,  the  viscosity  taking  into  account  fluid  particles  interactions  needs  to  be 
parameterized. 

The  governing  equation  of  kinetic  energy  for  subgrid-scale  flow  considering 
the  interaction  of  fluid  and  particles  can  be  derived  as  follow. 
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From  EQ  (5),  however,  the  third  order  fluctuating  terms  were  neglected,  Yuu 
(1997)  proposed  a two  way  coupling  SGS  model  applying  a local  equilibrium 
assumption  expressed  as  EQ  (6)  that  energy  production,  dissipation  and 
interaction  of  fluid-particles  are  in  balance. 
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aT.,  + 1 (7,  dx 
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From  EQ  (6),  the  viscosity  taking  into  account  of  fluid-particles  interaction  was 
derived  and  expressed  as  EQ  (7) 
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However,  it  is  difficult  to  decide  the  coefficient  CvT  dynamically  in  Yuu'  s 


model  because  the  filter  width  A are  mixed  in  EQ  (7).  In  this  study,  we  take  the 
place  of  EQ  (6)  Tu  with  EQ  (8)  and  get  EQ  (9)  as  follows 
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The  SGS  turbulence  energy  k can  be  modeled  as  follows  according  to 
dimensional  analysis 

*=P^(2V,)=^jA’|s|!  (10) 

Substitution  of  EQ  (10)  into  the  third  term  denominator  of  EQ  (9)  and  put  it  in 
order,  the  new  eddy  viscosity  is  then  given  by 
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Because  the  filter  width  A appears  at  only  one  place  in  the  EQ  (11)  similar  to  the 
standard  Smagrinsky  model  of  single-phase  turbulent  flows,  the  coefficient 
Cvr  of  proposed  model  can  be  easily  decided  by  Germano’s  (1991)  dynamic 
procedure. 

3.  A proposal  of  new  dynamic  SGS  model  of  Two  Way  Coupling 

Now  substitute  EQ  (11)  into  EQ  (2)  to  get  the  new  dynamic  eddy  viscosity 
subgrid-scale  stress  expressed  as  EQ  (12),  which  takes  into  account  of  the 
interaction  of  particles-fluid  SGS  components  in  LES. 
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According  to  the  Germano’s  dynamic  procedure  (l),  the  test  filter  level  SGS  stress 
can  be  defined  as 
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Where,  A is  the  characteristic  filter  width  associated  with  the  grid  level  filtering 
operator,  and  As  is  the  filter  width  associated  with  the  test  level  filtering 

operator.  Considering  the  similarity  between  the  SGS  stresses  at  the  grid  and  test 
levels,  which  are  modeled  using  the  same  functional  expression,  it  was  assumed 
that  the  coefficient  C\,T  or  CmS  in  EQ  (12)  and  EQ  (13)  are  same,  then  grid  filter 


stress  Ty,  test  filter  stress  T:j  and  resolved  turbulent  stress  can  be  related  in 
algebraic  relation  EQ  (14) 

Tj  -fj  = iji/y -ufij  = 2C^2 A[^^,upi, N)Si}  -«2^(«;,up,.,yV)^  (14) 

From  EQ  (14),  CvT(x,y,z,t)  can  be  obtained  in  principle.  The  quantity  in 
square  brackets,  however,  can  become  zero,  which  would  make  CvT 
indeterminate  or  ill-conditioned.  For  the  channel  flow,  therefore,  it  was  assumed 
that  CvT  is  only  a function  of  y and  t . The  average  of  both  sides  of  EQ  (14) 


are  taken  over  a plane  parallel  to  the  wall  to  yield 
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The  new  dynamic  eddy  viscosity  subgrid-scale  stress  model,  which  takes  into 
account  of  the  interaction  of  particles-fluid  SGS  components  in  LES,  is  then 
given  by 
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4.  Simulation  overview 
4.1  Calculation  of  fluid  phase 

The  particle-laden  turbulent  flows  between  plane  channels  driven  by  uniform 
pressure  gradient  and  particles  gravity  were  calculated  by  large  eddy  simulations 
at  Reynolds  numbers  based  on  friction  velocity  and  channel  half-width  of  644. 
Van  Driest  damping  function  model  and  dynamic  SGS  models  with  (the  present 
model)  or  without  (Germano’s  model)  couplings  of  SGS  components  were 
applied  respectively.  The  governing  equation  EQ  (1)  and  continuity  equation 
were  solved  numerically  by  SMAC  method  on  a staggered  grid.  Second-order 
central  difference  scheme  was  used  for  the  advection  and  diffusion  terms,  and 
second-order  Adams-Bashforth  method  was  adopted  for  time  advancement.  The 
Poisson  equation  for  pressure  was  solved  using  ICCG  method.  The  flow  was 
resolved  using  32  X 64  X 32  grid  points  in  the  x,  y and  z directions,  respectively. 
The  channel  domain  for  the  calculation  was  7r5X26X7r5/2.  For  fully 
developed  channel  flow,  periodic  boundary  conditions  for  the  dependent 
variables  were  applied  in  the  streamwise  and  spanwise  directions,  whereas  the 
no-slip  condition  was  applied  on  the  channel  walls. 
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4.2  Calculation  of  particle  phase 

The  motion  of  particles  was  integrated  using  second-order  Adams-Bashform  in 
time,  and  third  -order  Lagrange  polynomials  were  used  to  interpolate  the  fluid 
velocity  to  the  particle  position  since  it  is  only  by  chance  that  a particle  is  located 
at  a grid  point  where  the  Eulerian  velocity  is  available.  For  particles  that  moved 
out  of  the  channel  in  the  streamwise  or  spanwise  directions,  the  periodic 
boundary  conditions  were  used  to  introduce  them  into  the  computational  domain. 
The  channel  walls  were  perfectly  smooth  and  a particle  was  assumed  to  contact 
the  wall  when  its  center  was  one  radius  from  the  wall.  Elastic  collisions  were 
assumed  for  particles  contacting  the  wall. 

The  initial  condition  of  single  phase  Eulerian  velocity  field  was  given  by  a 
statistically  developed  solution.  Then  the  particles  were  assigned  to  random 
locations  throughout  the  cannel,  where  the  initial  particle  velocity  was  assumed 
to  be  the  same  as  the  fluid  velocity  at  the  particle  location.  Similar  to  the  fluid 
flow,  statistics  of  the  particle  velocity  were  averaged  over  two  homogeneous 
directions,  both  channel  halves  and  time.  To  distinguish  the  effect  of  two  way 
coupling,  the  70  p,  m Copper  particles  were  blended  at  mass  loading  ratio  1.0. 


5.  Calculation  results 


Figure  1 . Profile  of  streamwise  mean  velocity  of  fluid  Figure  2.  Turbulence  intensity  of  fluid  in  streamwise 

The  numerical  simulations  were  carried  out  using  Van  Driest  wall  function  and 
the  present  dynamic  SGS  model,  while  the  SGS  coupling  was  taken  into  account 
or  not  respectively.  The  streamwise  mean  velocities  of  fluid  are  shown  in  Figure 
1,  the  root-mean-square  velocity  fluctuations  of  fluid  in  streamwise  are  shown  in 
Figure  2.  As  may  be  observed  in  figures  above,  in  spite  of  the  SGS  coupling  is 
taken  into  account  or  not,  the  mean  velocities  are  predicted  larger,  and  the 
root-mean-square  velocity  fluctuations  are  also  slightly  larger  using  proposed 
dynamic  SGS  models  than  using  Van  Driest  wall  function.  Since  this 
characteristic  is  similar  to  single-phase  turbulent  flow  using  standard  dynamic 
SGS  model,  we  can  say  that  the  present  SGS  model  has  proper  asymptotic 
behavior  near  the  wall  without  the  use  of  ad  hoc  damping  function  in  case  of 
multiphase  turbulent  flow.  As  shown  in  Figure  2,  the  turbulence  intensities  of 
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particle-laden  flow  are  clearly  decreased  than  ones  of  particle-unladen  flow  for 
the  coupling  of  GS  components  of  particle-laden  flow,  but  the  turbulence 
intensities  of  particle-laden  flow  are  hardly  decreased  when  the  SGS  coupling 
was  taken  into  account. 


Figure  3.  Profile  of  mean  model  coefficienls  Figure  4.  Profile  of  mean  SGS  eddy  viscosity 

The  eddy  viscosity  model  coefficients  that  were  calculated  by  proposed 
dynamic  SGS  model  are  shown  in  Figure  3.  In  the  channel  buffer  region,  the 
dynamic  SGS  model  coefficient  of  particle-laden  flow  is  smaller  than  that  of 
particle-unladen  flow.  On  the  contrary,  in  the  channel  center  area,  the  former  is 
higher  than  the  latter.  Since  the  effects  of  two-way  coupling  of  GS  and  SGS 
components  correspond  the  profile  of  turbulence  intensity  of  fluid  in  wall  normal, 
the  coefficient  of  proposed  dynamic  SGS  model  obtains  the  correct  asymptotic 
behavior  in  the  near-wall  region.  In  logarithm  law  region(y+~100),  where 
concentration  of  the  local  particles  is  lowest,  the  proposed  model  obtains  the 
fittest  value  Cm~0.1  in  single-phase  channel  turbulent  flows.  It  reflects  the  local 
structure  of  small  eddies,  showing  that  the  capability  of  proposed  model  is 
validated. 

The  mean  distributions  of  eddy  viscosity  are  shown  in  Figure  4.  In  the 
near-wall  region  (y+<30),  as  the  effects  from  SGS  component  to  GS  component 
are  relatively  small,  the  distribution  of  model  coefficients  correspond  the  profile 
of  eddy  viscosity  as  they  are.  While  in  the  logarithm  law  region,  though  the 
model  coefficients  are  roughly  same,  the  effect  of  two-way  coupling  of  GS 
components  is  significant,  which  may  result  from  the  difference  of  strain  velocity 
S.  As  a result  of  taking  into  account  of  two-way  coupling  of  SGS  component,  the 
model  coefficients  become  bigger,  which  is  agreed  with  the  tendency  that 
dissipation  of  turbulence  is  enhanced  in  high  frequency. 

The  cascade  profile  of  GS  turbulence  energy  with  or  without  SGS  component 
coupling  is  shown  in  Figure  5,  the  energy  cascade  becomes  slightly  stronger 
owing  to  SGS  coupling  as  expected.  This  indicates  that  the  contribution  of  SGS 
component  to  GS  component  is  shown  as  dissipation  of  GS  component. 
According  to  this  simulation  result,  the  efficiency  of  presented  dynamic  two-way 
coupling  SGS  Model  is  verified.  Intriguingly,  contrary  to  the  conjecture  of 
Yamamoto*5 } et  al.  (1998),  the  spatial  spectra  of  fluctuation  energy  becomes  a 
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little  stronger  with  SGS  coupling  than  that  without  it  as  shown  in  Figure  6. 


Figure  5.  Cascade  Profile  of  GS  turbulence  energy 
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6.  Concluding  remarks 

(1)  A new  dynamic  SGS  model  has  been  proposed  for  two  way  coupling  LES 
of  particle-laden  turbulent  flow  based  on  Yuu'  s SGS  model,  in  which  tda 
couplings  of  fluid-particles  SGS  component  on  fluid  turbulence  modulation 
was  taken  into  account  and  at  the  same  time  the  coefficient  of  proposed 
SGS  model  can  be  dynamically  obtained  as  the  calculations  progress  rather 
than  input  a priori. 

(2)  As  a result  of  taking  into  account  of  two-way  coupling  of  SGS  component, 
the  dynamic  model  coefficients  become  bigger,  and  in  the  channel  buffer 
region,  the  proposed  dynamic  model  coefficient  of  particle-laden  flow  is 
smaller  than  that  of  particle-unladen  flow,  on  the  contrary,  the  former  is 
higher  than  latter  in  the  channel  center  area. 

(3)  Numerical  simulation  results  using  the  new  dynamic  model  proposed 
reflected  the  local  structure  of  small  eddies  of  particle-laden,  showing  that 
the  capability  of  proposed  model  was  validated. 
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